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It has been shown recently that non-Abelian states and the spin-polarized and unpolarized ver- 
sions of the Abelian 331 state may have identical signatures in Fabry-Perot interferometry in the 
quantum Hall effect at filling factor 5/2. We calculate the Fano factor for the shot noise in a Mach- 
Zehnder interferometer in the 331 states and demonstrate that it differs from the Fano factor in the 
proposed non-Abelian states. The Fano factor depends periodically on the magnetic flux through 
the interferometer. Its maximal value is 2 x 1.4e for the 331 states with a symmetry between two 
flavors of quasiparticles. In the absence of such symmetry the Fano factor can reach 2 x 2.3e. On 
the other hand, for the Pfaffian and anti-Pfaffian states the maximal Fano factor is 2 x 3.2e. The 
period of the flux dependence of the Fano factor is one flux quantum. If only quasiparticles of one 
flavor can tunnel through the interferometer then the period drops to one half of the flux quantum. 
We also discuss transport signatures of a general Halperin state with the filling factor 2 -I- fc/(A; + 2). 
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I. INTRODUCTION 

Two-dimensional electron systems in a strong mag- 
netic field exhibit much beautiful physics and form many 
states of matter including numerous fractional quantum 
Hall liquids^. Some of their properties such as frac- 
tional charges of elementary excitations are well under- 
stood. Less is known about the statistics of quantum 
Hall quasiparticles. We know that gauge invariance re- 
quires fractionally charged particles to be anyons. At the 
same time, a direct experimental observation of anyonic 
statistics poses a major challenge. This challenge has 
recently attracted much attention because of a possibil- 
ity of non-Abelian anyonic statistics at some quantum 
Hall filling factors^. In contrast to "tamer" Abelian par- 
ticles, non-Abelian anyons change their quantum state 
after one particle makes slowly a full circle around other 
anyons. This property can be used for topological quan- 
tum computatioii^. Possible application as well as in- 
trinsic interest of such unusual particles have stimulated 
attempts to find non-Abelian anyons in nature. In partic- 
ular, a possibility of non-Abelian statistics was predicted 
at filling factors 5/2 and 1 jl^—. However, the nature of 
the quantum Hall states at those filling factors remains an 
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FIG. 1: Schematic picture of an anyonic Mach-Zehnder inter- 
ferometer. Arrows indicate propagation directions of the edge 
modes on Edge 1 (from source S\ to drain D\) and Edge 2 
(from source 5*2 to drain Di). Quasiparticles can tunnel be- 
tween the two edges through two quantum point contacts, 
QPCl and QPC2. 



open question with theoretical proposals including both 
Abelian and non-Abelian states^"—. 

Numerical simulationsi^"— with small systems of 8- 
20 electrons provide support to non-Abelian Pfaffian 
and anti-Pffafian states. This support is however not 
unanimous, see, e.g., Ref. [igI At the same time, re- 
cent experimentsii"— suggest an unpolarized state at 
V = 5/2. Zero spin polarization is incompatible with 
either Pfaffian or anti-Pfaffian states. Some of the 
experiments can be understood in terms of disorder- 
generated skyrmions^S in a Pfaffian state. However, such 
an explanation does not apply to the most recent opti- 
cal experimenti^ and the simplest interpretation of the 
existing limited experimental data is in terms of zero 
polarizatio n^ ^1^^ . The simplest unpolarized state is the 
Halperin 331 stated. Note that the results of the tunnel- 
ing experiment^ are compatible with both 331 and anti- 
Pfaffian states. In a closely related problem of quantum 
Hall bilayers at filling factor 5/2, numerics supports the 
existence of both 331 and Pfaffian states, separated by a 
phase transitions^. Thus, it is important to find a way to 
identify and distinguish from each other the 331, Pfaffian 
and anti-Pfaffian states. 

The key difference lies in non-Abelian quasiparticle 
statistics of the Pfaffian and anti-Pfaffian states versus 
Abelian statistics in the 331 state. In order to determine 
the statistics of anyons at v = 5/2 and 7/2 several ex- 
periments were proposed and some of them have been 
or are being currently implemented (for a review, see 
Ref. [25I ). Despite those efforts the statistics in the 5/2- 
and 7/2-states remains an open question. One of the 
issues concerns ambiguities in the interpretation of the 
experimental data. In particular, the most elegant and 
conceptually simple approach to detecting non-Abelian 
anyons is based on Fabry-Pcrot interferometry^^"— . It 
was found recently^"— that a Fabry-Perot interferom- 
eter may produce identical interference and Coulomb 
blockade patterns in transport experiments with the non- 
Abelian Pfaffian and anti-Pfaffian states and the spin- 
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polarized^iSrii and unpolarized^ versions of the Abelian 
331 state. Similarity between the patterns is only present 
in the case of the exact or approximate symmetry be- 
tween two quasiparticle flavors in the 331 states and one 
may hope to remove such symmetry by some perturba- 
tion or a change in the conditions of the experiment^. 
However, in the absence of an estabhshed theory of 5/2- 
and 7/2-states, it is hard to tell whether the flavor sym- 
metry is or is not present and if a particular perturbation 
would make it possible to distinguish Abelian and non- 
Abelian states. Thus, it is desirable to have another ap- 
proach to interferometry such that Abelian states would 
never mimic non- Abelian. 

In this paper we show that in an anyonic Mach- 
Zehndcr interferometer—"—, signatures arc different for 
the Pfafhan and anti-PfafSan states on the one hand and 
the polarized and unpolarizcd 331 states on the other 
hand. This conclusion is true both in the presence and 
absence of the flavor symmetry. We calculate the low- 
temperature zero-frequency noise in the interferometer 
in the limit of weak tunneling through the device. In 
such case the noise is related to the current as 5 = 2e*I, 
where 2e* is the Fano factor. The Fano factor exhibits 
a periodic dependence on the magnetic flux. Its maxi- 
mal value as a function of the flux was calculated for the 
Pfaffian and anti-PfafHan states in Ref. [13 and equals 
2e* = 2 X 3.2e, where e is an electron charge. We show 
that in the 331 states with flavor symmetry the maxi- 
mal value of the Fano factor is 2e* = 2 x 1.4e. In the 
absence of the symmetry the Fano factor can reach the 
maximal value of 2 x 2.3e in the 331 states. Thus, the 
maximal Fano factor gives an unambiguous way to dis- 
tinguish the Abelian 331 quantum Hall liquids from the 
proposed non- Abelian states. We also calculate the elec- 
tric current through the interferometer as a function of 
the flux and voltage. The results for the 331 state dif- 
fer from the case of the Pfaffian state but the difference 
between the two I-V curves is small. 

The paper is organized as follows. First, we briefly 
discuss the 331 states. Next, we review the structure 
of the anyonic Mach-Zehnder interferometer in Sec. IIIII 
We calculate the zero-temperature current in Sec. II VI and 
zero-temperature zero- frequency shot noise in Sec.|Vl Wc 
summarize our results in Sec. lVIl The Appendix contains 
a detailed discussion of general Halperin states at flUing 
factor ly = 2 + k/{k + 2). 



II. STATISTICS IN THE 331 STATE 

A review of different proposed states for the filling fac- 
tor 5/2 can be found in Ref. 0- Here wc summarize the 
properties of the 331 states^"—. 

Two different variants of the 331 state are described 
in the literature. The spin-unpolarized version was in- 
troduced in Ref. [^. It can be understood as a bilayer 
state with two spin components playing the role of the 
layers. The filling factor in each layer is 1/4. The spin- 



polarized version^"— can be described in terms of the 
condensation of the charge-2e/3 quasiparticles on top of 
the Laughlin = 1/3 state. The two states differ in 
many respects but have the same key features: the topo- 
logical order and the statistics of quasiparticles^. Since 
anyonic interferometry is only sensitive to the quasipar- 
ticle statistics, it cannot distinguish the two states. Be- 
low we describe the statistics with the help of the K- 
matrix formalism, Ref. H^. We only focus on the half- 
filled Landau level. Integer edge channels formed in the 
lower completely filled Landau levels are unimportant for 
our problem since the transport through the interferome- 
ter is dominated by the tunneling of fractionally charged 
excitations on top of the half-filled Landau level. 

The JC-matrix formalism encodes the information 
about quasiparticles in terms of a matrix K and a charge 
vector t. Each elementary excitation is described by a 
vector 1„ with integer components. In the 331 states all 
edge modes propagate in the same direction. In such case 
the scaling dimensions of quasiparticle creation and anni- 
hilation operators are independent of the interactions be- 
tween the modes and are given by /i„ = l„i^~^l^. We are 
only interested in the most relevant quasiparticle opera- 
tors. The quasiparticle charge Q = — et„iC~^lJ. Trans- 
port through the Mach-Zehnder interferometer depends 
on statistical phases accumulated by the wave function 
when one particle makes a full circle around another. The 
phase accumulated by particle 1 moving around particle 
2 equals 612 = 27rli7^-il|'. 

The spin-unpolarized 331 stated can be described as a 
bilayer state with ~ i in each layer. The X-matrix is 

K^Cli). ,1, 

The charge vector t=(l,l). The two most relevant quasi- 
particles are characterized by the Z-vcctors li = (1, 0) and 
I2 = (0, 1). Both particles have charge e/4. This elemen- 
tary excitation charge agrees with experiments^ii^Siii. 
When a particle makes a circle around an identical 
particle it accumulates the phase 0ii = ^22 ~ 37r/4. 
The mutual statistical phase of two different particles 
is 012 ~ — 7r/4. We will also need to know what 
phases are accumulated when a charge-e/4 quasiparticle 
qo moves around a composite anyon built from several 
e/4-quasiparticles qi,...,qk- In the Abelian 331 state 
such statistical phase is simply the sum of mutual statis- 
tical phases of go with each of the particles. 

The spin-polarized 331 stated is formed by the con- 
densation of the charge-2e/3 quasiparticles on top of the 
Laughlin = 1/3 state. This state is characterized by 
the A'-matrix 

'<-{-.-:)■ 

and the charge vector t = (1,0). Calculations of the 
electric charges and statistical phases are the same as 
above. The two most relevant elementary excitations 
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carry charges e/4 again and arc characterized by vec- 
tors li = (0,1) and h = (1,-1). All statistical phases 
are the same as in the spin-unpolarized 331 state. 

In fact, the two versions of the 331 state are topologi- 
cally equivalenlj^, since the two iiT-matrices satisfy 

and the charge vectors satisfy 

(1,1) = (1,0)1^^ (4) 

This explains why the two states have the same quasi- 
particle charges and braiding statistics. 

III. MACH-ZEHNDER INTERFEROMETER 

Figure [T] shows a sketch of the Mach-Zchndcr geometry. 
Because of the bulk energy gap, the low-energy physics 
is determined by chiral edge modes. Charge flows along 
Edge 1 from source 5i to drain Di and along Edge 2 
from source 52 to drain D2- Quasiparticles tunnel be- 
tween the two edges at quantum point contacts QPCl 
and QPC2. If one keeps at a positive voltage V and 
the other source and drains are grounded then there is 
a net quasiparticle flow into Edge 2 and a net tunneling 
current. The current is measured at drain D2. 

As discussed above, there are two flavors of charge- 
e/4 quasiparticles in the 331 state. Let us denote their 
topological charges (or flavors) as a and b. Since the 
most relevant quasiparticle operators create particles of 
these two types, we consider only the tunneling of e/4- 
quasiparticlcs with flavors a and b below. We focus on the 
limit of small tunneling amplitudes between the edges. In 
such case the problem can be accessed with perturbation 
theory. We denote the small tunneling amplitudes at the 
two point contacts as F^, where /c = 1, 2 is the number of 
the point contact and x = a,b is the topological charge 
of the tunneling quasiparticle. The tunneling rate from 
Edge 1 to Edge 2 can be found from the Fermi golden 
rule. It depends on the tunneling amplitudes and on the 
phase difference for the quasiparticles which follow from 
SI to D2 through QPCl and QPC2. The latter consists 
of two contributions: the Aharonov-Bohm phase due to 
the external magnetic field and the statistical phase ac- 
cumulated by a quasiparticle making a full circle around 
the "hole" in the interferometer. The statistical phase is 
determined by the total topological charge that tunneled 
previously between the edges. Indeed, the total topolog- 
ical charge of Edge 2 can only change during tunneling 
events at QPCl and QPC2. Charge exchange with the 
Fermi-liquid drain D2 and source S2 cannot affect the 
topological charge of the edge. Taking into account that 
the edges are chiral we see that all topological charge 
that previously tunneled into Edge 2 accumulates inside 
the loop QPC1-A-QPC2-B-QPC1. Wc will denote that 
loop as £ below. Certainly, the accumulated topological 



charge only assumes a discrete set of values and hence 
changes quasiperiodically as a function of time. We find 
the following tunneling rate from edge 1 to edge 2 for a 
quasiparticle of flavor x (cf. Ref. [37| ) 

<d - ridFf |2+|Ff n+(r2F^rre"^»-+'*==^+c.c.), (5) 

where ri{V,T,x) and r2{V,T,x) depend on the voltage, 
temperature and quasiparticle flavor, d is the topologi- 
cal charge trapped in the interferometer before the tun- 
neling event, (j)xd the statistical phase discussed in the 
previous section, and the Aharonov-Bohm phase 0mag = 
7r$/(2$o) is expressed in terms of the magnetic flux $ 
through the loop £ and the flux quantum $o = he/ e. ri_2 
cannot be calculated without a detailed understanding of 
the edge physics. Fortunately, we will not need such a 
calculation to determine the main features of the current 
and noise in the interferometer. Eq. ([5]) assumes that the 
two tunneling amplitudes are small and hence the mean 
time between two consecutive tunneling events is much 
longer than the time spent by a tunneling quasiparticle 
between the point contacts. At a nonzero temperature, 
quasiparticles are allowed to tunnel from Edge 2 to Edge 
1 which has a higher potential. The corresponding tun- 
neling rate w~ — cxp[~eV/{AkBT)]w^ ^ is connected 
to Eq. ([5]), by the detailed balance principle. Here 
d+x is the topological charge of the fusion of anyons with 
topological charges d and x, i.e., the topological charge 
of Edge 2 before the tunneling event. In what follows we 
concentrate on the limit of low temperatures and neglect 
w~ . 

Our main focus will be on the situation with the ex- 
act or approximate flavor symmetry, i.e., we assume that 
the tunneling amplitudes F^ and coefficients rk{V,T,x) 
in (O do not depend on the flavor x. Since F's and 
r's only enter the transition rates in the combinations 
ri{V, T, a;)|Ff_2p and r2{V, T, a;)Ff F|*, the results do not 
change if the tunneling amplitudes depend on the flavor 
but the above combination are the same for x = a and 
b. As discussed in the introduction, Fabry-Pcrot intcr- 
ferometry cannot distinguish the 331 states with flavor 
symmetry from non-Abelian states. We will see below 
that no such issue exists for Mach-Zchnder interfcrom- 
ctry. Wc will also briefly discuss signatures of the 331 
states without flavor symmetry in the shot noise in a 
Mach-Zchnder interferometer. 



IV. ELECTRIC CURRENT 

We arc now in the position to calculate the tunnel- 
ing current /. From now on, we will only use the lan- 
guage of the spin-unpolarized 331 state in which flavor 
a can be understood as the Z-vector a = (1,0) and b as 
b = (0,1). From Sec. |TT] we know how to evaluate the 
statistical phase in The tunneling rate depends on 
the topological charge d ~ (m, n) trapped in the inter- 
ferometer. Each tunneling event of a particle with flavor 
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(a) 




FIG. 2: Possible states of a Mach-Zehnder interferometer in the 331 state. Panel (a) shows a general case, eight possible states 
labeled by topological charges and the transition rates between them. Arrows show the allowed transitions at zero temperature. 
Solid blue lines represent tunneling events involving quasiparticles of flavor a, and dashed black lines represent tunneling events 
involving particles of flavor b. Special cases with — p\ = pk and = are illustrated in Panels (b) and (c) respectively. 



a changes d — (m+ 1, n). If a particle of flavor b tunnels 
then d — >■ (m, n + 1). 

Any anyon has a trivial mutual statistical phase 27rfc 
with an electron^. Combining this condition vifith the 
knowledge of the electron charge one finds that the l- 
vectors (1,3) and (3,1) describe electrons. Thus, the 
topological charges d and d' = d + ni(l,3) + n2(3,l) 
can be viewed as identical since anyons accumulate iden- 
tical topological phases moving around charges d and d'. 
One can easily see that in a 331 state, the trapped topo- 
logical charge falls into one of eight equivalence classes 
which mark eight possible topological states of the area 
enclosed by the loop £,. Fig. [2{a) shows the 8 states and 
possible transitions between them due to anyon tunneling 
between the edges at zero temperature. The transition 
rates shown in Fig. [5] are given by the equation 

pI = A'^ll + u"" cos(7r$/(2$o) + T^k/A + 5"=)] (6) 

where x = a ov b, and k = 0, 1, . . . , 7. The parameters 
A- = ridFf |2-H|rf = 2|r2FfFi|/[ri(|Ff |2-H|Ff |2)] 

and (5^ = arg(r2Ff Fg*). In the presence of the flavor 
symmetry the diagram in Fig.[2ja) simplifies to Fig. [2Kb) 
with six instead of eight vcrtexes. Two pairs of vertexes 
in Fig. [2] (a) are merged into single vertexes in Fig.[2](b). 
This is legitimate due to numerous equalities between 
tunneling rates. In Fig. [2] (b), we use the notation pk = 
p'^ =p\. We will also denote A A°- ^ A'' , u = u°- = u'' 
and 6 = 6" = 6^ in the presence of the flavor symmetry. 
Fig.^c) illustrates another simple limit in which = 0. 




FIG. 3: Flux-dependence of the tunneling current in the 331 
and Pfafhan states. We set = 1, = m'' = 1, and 
_ gb _ Q £|_^j. curves for the 331 state. Different curves 
correspond to different values of 7 = A'' /A"' in the 331 state. 
The curve for the Pfaf&an state is plotted according to Eq. (8) 
in Ref. [3^ with rf^ = rf^ = 1 and Fi — F2 such that the max- 
imum matches the maximum of the curve for the 331 state 
with 7 = 1. 



The transition rates can be used to write down a ki- 
netic equation for the probabilities fd of different topo- 
logical charges d trapped in the interferometer. In terms 
of the distribution function the average tunneling cur- 
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rent between the edges 

d 

where d goes over the eight possible states in Fig. ^a). 
The distribution function satisfies the steady state equa- 
tion 

^^^^ J2 ifd-^'^td-x + fd+xW-^^J 

X—a.b 
x—a^b 



and the normahzation condition J2d fd = ^- In Eq. ([8|), 
{d + x) means the topological charge of the fusion of 
anyons with topological charges d and x. The fusion of x 
and d — X has topological charge d. In the absence of the 
flavor symmetry, the tunneling current can be found with 
a lengthy but straightforward calculation from Eqs. (O 
and ([S]). In the presence of the flavor symmetry, the low- 
tcmpcrature current can be easily calculated from Eqs. 
(|7l8p in the picture with six distinct topological charges, 
Fig. [21(b). One obtains the following result 



eA 



+ ^{1 - cos(27r$/$o + 4(5)) 



4V2' 
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(1 - ^)(1 - cos(27r$/$o + 4(5)) - ^ sin(27r$/$o + 4(5) 
I 



(9) 



This formula is quite similar to the expression for the 
current in the Pfaffian state (cf. Eq. (8) in Ref. l36l ). 
The similarity originates from the similarity of the dia- 
gram Figl^fb) with a corresponding diagram in the Pfaf- 
fian stated. They have the same topology with six ver- 
texes including two "cross-roads" . Still, in contrast to 
the Fabry-Perot case, the expressions for the current are 
not identical for the 331 and Pfaffian states. Similar to 
the Pfaffian stated, the expression for the current for 
the opposite voltage sign can be obtained by changing 
both the overall sign of the current and the sign be- 
fore sin(27r$/$o + 4(5) in the denominator. Thus, the 
I-V curve is asymmetric just like in the Pfaffian case. 
The current depends periodically on the magnetic flux 
with the period $0 in accordance with the Bycrs-Yang 
theorem^. We do not include an analytical expression 
for the current in the absence of the flavor symmetry 
since it is lengthy. /(<&) is plotted in Fig. 3 for the Pfaf- 
fian and 331 states for different values of 7 = A'^/A'^ at 
u = 1 which maximizes the visibility of the Aharonov- 
Bohm oscillations. 

Another simple limit corresponds to the situation with 
= (Fig. IHc)). In that case, /($) has a reduced 
period $0/2. The period reduction can be understood 
from the diagram Fig.[2jc). The system can return to its 
initial state only after eight tunneling events instead of 
four in Fig. [2] (b). This means a transfer of the charge 
2e in each cycle of tunneling events. Such "Cooper pair" 
charge agrees with the "superconductor" periodicity. 

V. SHOT NOISE 

Shot noise measurements have helped to determine the 
charges of elementary excitations at several quantum Hall 



filling factors including 5/3^i^i^. In the Mach-Zehnder 
geometry, zero-frequency noise contains also information 
about the quasiparticle statistics. Below we calculate 
the noise in the 331 states in the limit of weak tunneling. 
We assume that the temperature is much lower than the 
voltage bias. In such case it is possible to neglect the 
temperature and perform calculations at T = 0. 

Shot noise is defined as the Fourier transform of the 
current-current correlation function, 

/+00 
(/(0)/(i) + /(i)^(O)) exp{iu}t)dt. (10) 
-00 

Below we consider w = only. In the weak tunneling 
limit, the noise can be expressed as S* = 2e*/, where 
2e* is known as the Fano factor, e* can be understood 
as an effective charge tunneling through the interferom- 
eter. We will see that the Fano factor is different in the 
Pfaffian state and the 331 states with or without the fla- 
vor symmetry. In the Mach-Zehnder interferometer, the 
Fano factor 2e*($) exhibits oscillations as a function of 
the magnetic flux. For Laughlin states, the maximal e* 
can never exceed l.Oe, while in the Pfaffian state e* can 
be as large as 3.2e, Ref. [s^- We will see that the maximal 
Fano factor is lower in the 331 state than in the Pfaffian 
state. 

A zero-frequency shot noise can be conveniently con- 
nected with the fluctuation of the charge transmit- 
ted through the interferometer over a long time t (cf. 
Ref. [13): 

S/2 = lim {6Q^{t))/t, (11) 

where 5Q{t) is the fluctuation of the charge Q{t) that has 
tunneled into the loop £ during the time t. The average 
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FIG. 4: The msiximal e* as a function of 7 and 5. 



current / = \irat^ao{Q{t)) /t- We will use the genera ting 
function method developed in Ref. [13 (see also Ref. l38h 
to calculate the zero-frequency shot noise and the Fano 
factor. Suppose for certainty that the initial state inside 
the loop £ has topological charge (0,0). Our results will 
not depend on this choice of the initial topological charge. 
The electric charge that had tunneled through the inter- 
ferometer, Q{t), is zero at t = 0, (5(0) = 0. The system 
evolves over a time period t according to a kinetic equa- 
tion with the transition rates from the diagram Fig.[2ja). 
The probability of finding the interferometer in the state 
with the topological charge d on Edge 2 and the transmit- 
ted electric charge Q{t) = fc| will be denoted as fd,k{t) 
below. For a given d, only certain fc's are possible (e.g., 
if d = (0, 1) then k can only be 4ri -I- 1 with an integer n). 
The distribution function satisfies the kinetic equation 

dfd,k{t) t f + , r - X 

x—a.b 

- E (12) 

X—a.b 



Let US now define the generating function 



/d(s,0 = E/'^.'^W^'' 



(13) 



where the summation extends over all possible k's for a 
given d. One can express the average transmitted charge 
{Q{t)) — eJ2k kld,k{t)/^ and its fluctuation as 



(14) 



and 



w^"»^(i)'(l'|i:«».')) 



where the angular brackets stay for the average with re- 
spect to the distribution function fd,k{t)- It is easy to 
see that fd{^,t) equals the probability fd{t) introduced 
in Sec. IIVI The kinetic equation then reduces to 



x—a,b 

- E fdiw^A + W~,d) 
x—a,b 



(16) 



The above equation can be rewritten in a matrix form, 

f{s,t) = M(s)/(s,i), where M{s) is a time-independent 
8x8 matrix. We can solve the above linear equation and 
obtain 



fd{s,t) = Y,uds)^^''-'^'m), 



(17) 



where Xi{s) arc eigenvalues of M{s), £^d.i constants and 
Pi{t) polynomials, Pi = 1 for nondcgenerate eigenvalues 
\i . The solution can be further simplified with the help of 
the Rohrbach theorem^. Indeed, from (fTH]) we see that 
M (s) has negative diagonal elements, non-negative off- 
diagonal elements, and the sum of the elements of each 
column is zero at s = 1. In such situation the Rohrbach 
theorem applies. According to the theorem, the eigen- 
value of M(s = 1) with the maximal real part is non- 
degenerate and equals zero. All other eigenvalues have 
negative real parts. For s close to 1, the maximal eigen- 
value Ao(s) should also be nondcgenerate by continuity 
of the eigenvalues as functions of s. Thus, the Ao(s) term 
in Eq. (flT)) dominates at large time t and /o(i) = 1- 

According to p4)l . (fT5|) and (flT)) together with the con- 
servation of probability fd{l,t) = 1, we have 



^ = ^a;,(i) 



(e/4)2 [A;,'(1) + A[,(1)] 
^[1 + A^,'(1)/A^,(1)] 



(18) 



The derivatives of Aq (s) can be evaluated in the following 
way. Define a function 



G(s, A) = det(M(s) - \E) 



(19) 



(15) 



where det is the determinant of a matrix and E is the 
identity matrix. We know that Ao(s) satisfies the equa- 
tion G(.s, Ao(s)) = 0. Then by differentiating this equa- 
tion one gets 

'^o(l) = -Gs/Ga|(s,a)=(i,o) 

-^o(l) ^ —{GssG\ + G\\G\ — 2GsaGsGa)/G'^|(s_a)=(i.o)i 

(20) 

where Gs and Ga are the first derivatives of G(s, A) with 
respect to s and A, and G^g, Gsa and Gaa are the three 
second derivatives of G(s, A). Then we can easily evaluate 
the shot noise and the Fano factor. 
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When the symmetry between the two flavors of quasi- 
particles exists the expression for the Fane factor at zero 
temperature simphfies: 



A;,(1) = 16^/(4 



P2fe+1 . 



P2k+i 



16 
1 



A:=0 

-y 



4A 



P2fc+1 
P2fc+4 



1 P2fc+1 



16A2 ^ P2k+4 



4 P6 P2 J34 Po' 



(21) 



where the identity +Pk+4: ~ '2 A is used to simplify the 
expressions. Like the current, the Fano factor e* is a pe- 
riodic function of the magnetic flux $ with the period $o- 
Numerical analysis of Eq. (j2ip shows that maximal e* is 
1.4e, greater than l.Oe for Laughlin states but smaller 
than 3.2e for the Pfaffian state. The maximal Fano fac- 
tor is achieved at u = 1. This value of u corresponds 
to Fi = F2 and ri = r2. As can be seen from a renor- 
malization group picture, the latter relation between is 
satisfied at low voltages and temperatures, eV, T < hv/a, 
where v is the velocity of the slowest edge mode and a 
the distance between the point contacts along the edges 
of the interferometer. The values of F^ can be controlled 
with gate voltages. Note that m = 1 is the maximal pos- 
sible value of u. Indeed, u > 1 would result in negative 
probabilities (l6|) at some values of the magnetic flux. 

An interesting situation emerges in the special limit 
when only particles of one flavor a or 6 can tunnel (i.e., 
A° = or = 0). The interferometer can be tuned 
to that limit with the following approach: one keeps the 
filling factor ly = 5/2 in the gray (green online) region 
in Fig. 1 and ly — in the white region. In addition, a 
narrow region with a filling factor 7/3 is created along the 
edges. In such situation, quasiparticles tunnel through 
the 5/2-liquid between the interfaces of the v = 5/2 and 
1^ — 7/3 regions at the tunneling contacts. Since the 
interface contains only one edge mode, only one type of 
quasiparticles can tunnel (for a detailed discussion of the 
interface mode sec Appendix, subsection 1). 

We will assume for certainty that = 0. This case 
is illustrated in Fig. [2] (c). The calculation of the zero- 
temperature current and Fano factor greatly simplifies in 
that limit. One finds expressions resembling the results 
for the Laughlin states^: 



2e 



2e 



[E^Pk? 



(22) 



(23) 



Combining the above equations with Eq. ^ one can see 
that in the special case of only one type of quasiparti- 
cles allowed to tunnel, the current and noise are periodic 



functions of the magnetic flux with the "superconduct- 
ing" period $o/2. This is certainly compatible with the 
Bayers- Yang theoreroH. The numerator of the fraction 
in Eq. (|23p is always smaller or equal to the denomina- 
tor. They become equal when one of the probabilities pk 
approaches zero. In that case the Fano factor is maximal 
and e* = 2e. As is clear from Eq. this can happen 
only at u = 1. Note that the period and maximal Fano 
factor are also $0/2 and 2 x 2e in the K = 8 Abelian 
stated. 

Finally, let us discuss the most general case with 
nonzero A" and A'' and no flavor symmetry. Generally, 
when the flavor symmetry is absent, e*{^) is a periodic 
function with the period $0- Our numerical results show 
that the maximal e* is achieved at = u'' = 1 for any 
choice of A"", A'', i5° and 5''. We calculated the maximal 
e* = e*jj3^^(7, S) as a function of the tunneling amplitude 
ratio 7 = j4^/yl° and phase difference S = 5^ — (5°. We 
find that the maximal value of e^^^ = 2.3e is achieved 
at 7 = 0.08 and S = —0.03. The dependence of e'^^^ on 
7 and 6 is shown in Fig. [H We see that at 7 = 0, e*^^^ 
is 2e, increases to 2.3e at small positive 7 and drops to 
1.4e at 7 = 1. 

In a general case, the current and noise depend on 
several parameters and the system may not be tuned to 
the regime with the maximal Fano factor 2 x 2.3e. At the 
same time, such tuning is possible in the flavor-symmetric 
case, in the case when only one flavor tunnels and in the 
Pfaffian state. In all those cases one just needs to achieve 
u — I which corresponds to Fi = F2. Making Fi_2 equal 
is straightforward: one just has to make sure that the 
current is the same when only QPCl or only QPC2 is 
open. Thus, in the absence of the flavor symmetry, the 
identification of the 331 states simplifies by operating the 
interferometer in the regime in which only one quasipar- 
ticle flavor can tunnel. 



VI. SUMMARY 

We have calculated the current and noise in the 
Mach-Zehnder interferometer in the 331 state and 
compared the results with those for the Pfaffian state. 
Note that the transport behavior is essentially the same 
in the Pfaffian and anti-Pfaffian states. The current 
dependence on the magnetic fiux turns out to be quite 
similar for the Pfaffian and 331 states. In both states 
the I-V curves are asymmetric. The states can be 
unambiguously distinguished with a shot noise measure- 
ment. In the Pfaffian state the maximal Fano factor is 
2 X 3.2e. In the 331 state the Fano factor cannot exceed 
2 X 2.3e. If the flavor symmetry is present than the 
maximal Fano factor drops to 2 x 1.4e. The difference 
of the predicted maximal Fano factors well exceeds 
the current experimental accuracy of 15 percent^. An 
interesting situation emerges, if quasiparticles of only 
one ffavor can tunnel through the interferometer. In 
that case the current and noise arc periodic functions 
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of the magnetic flux with the period 3>o/2. In a general 
case the period is $o- We show that it is possible to 
tune the interferometer to the regime with the period 
$0/2 in the 331 state. 
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Appendix A: Multi-component Halperin states 

In this appendix we consider a general multi- 
component Halperin stated with the filling factor v = 
2 + k/{k + 2). In the first section of the appendix we de- 
scribe the states. In the second section we calculate the 
current and noise in the situation in which quasiparticles 
of one flavor dominate transport through the interferom- 
eter. This is the case of main interest for k > 2. In 
the third section we consider a general situation. Below 
we ignore the lowest filled Landau level and concentrate 
on the fractional quantum Hall effect in the second Lan- 
dau level. Indeed, the most relevant tunneling operators 
involve only the fractional edge modes. 

1. Quasiparticle statistics and edge modes 

The multi-component Halperin stato^. with the filling 
factor ly = 2 + k/{k + 2) can be described by a fc x /c ma- 
trix (-R'muit)ij = 1 + 26ij. All components of the charge 
vector t equal 1 and the most relevant quasiparticles are 
described by vectors 1^ = (0, . . . , 0, 1, 0, . . . , 0), where the 
number 1 stays in position i. The charge of the elemen- 
tary excitations is q = e/{k + 2). The state of the inter- 
ferometer is described by the numbers ni, . . . , njt of the 
trapped quasiparticles of each of the k types. The statis- 
tical phase, accumulated by a particle of type Im going 
around the hole in the interferometer, is 

9 = 27r^ if^^Tip = 7r[7i,„ - n/{k + 2)], (Al) 
p 

where n = ^n^. Hence, the tunneling probability for a 
particle of type m is 

p = A™ l-fu„ cos( +d„,-f7rn^--— — ) , (A2) 
(K-t-2)<l'o k + 2 

where Am, Um and 6m are real constants. 

An alternative description of the same topological state 
of matter can be formulated in terms of single-component 
hierarchical states. The starting point is the = 1/3 
Laughlin state. Condensation of quasiparticle pairs on 
top of the Laughlin state gives rise to the v = 2/(2 -I- 2) 



hierarchical state. Condensation of its quasiparticle pairs 
results in the ~ 3/(3-1-2) state and so on. The ap- 
propriate X-matrix has the following nonzero ele- 
ments: Kii = 3, Knn = 4 (/c > n > 1), if„,„_i = 
Kn-i n = ^2. The /C-matirx ii'muit expresses via ifh as 
i^muit = WK^W^, where W^, = 1, if z + j < (fc + 1), 
and Wij = 0, if i + j > k + 1. The charge vector 
reads t = (1,0, ...,0). The /-vectors of the most rel- 
evant quasiparticle excitations are 1 = (1, — 1, 0, . . . , 0), 
(0,1,-1,0,...,0),..., (0,...,0,1,-1), (0,...,0,1). 

In the main part of the article we considered k ~ 2 and 
focused on the flavor-symmetric case. This was justified 
for k = 2 for two reasons. First, in the flavor-symmetric 
case, the Fabry-Perot interferometry cannot distinguish 
the Pfaffian and 331 states. Second, for the unpolarized 
331 state the role of the flavors is played by the electron 
spin and an approximate symmetry between two spin 
projections in quantum Hall systems may give rise to the 
flavor symmetry. At fc 7^ 2 neither reason applies. In- 
deed, the Fabry-Perot interferometry can distinguish the 
states with fc ^ 2 from proposed non-Abclian states^ 
and k > 2 different flavors cannot be reduced to different 
spin projections. In the absence of the flavor symmetry 
one expects different tunneling probabilities for differ- 
ent quasiparticle types. As explained below, we expect, 
in general, that for one quasiparticle type the tunneling 
probability is much higher than for the other quasiparti- 
cles. 

To understand why this happens we need to identify 
edge channels in a ^ = fc'/(fc + 2) quantum Hall liquid. 
For this purpose we need to diagonalize the first term 
in the edge action^, Lk = J dtdx[K\i^ijdt(t>idx4>j — 
Vijdx4>idx4'j]- Note that we use the language of hier- 
arhical states. The diagonalization can be accomplished 
with the new variables $„ = (71 -I- 2)(/)„ — {n + l)</)„+i, 
n = 1, . . . , fc — 1 and $fe — {k + 2)(j)k- The action assumes 
the form 

Lk^Y I ^^^^tll A;°5t$„9,$„ - %^x^^^x^J], 
^ n=l 

(A3) 

where K° = 2/[{n + l){n + 2)]. The operator V« = 
exp(«4>„) creates an excitation with one electron charge. 
The coefficients and the charge created by V'n are 
independent of k. This allows for an easy description 
of the interface between the v = k/{k + 2) and = 
(fc— l)/(fc-|-l) liquids. Indeed, the action for the interface 
assumes the following form 

n=l Ti = l 

-|- electrostatic interaction -|- interchannel tunneling, 

(A4) 

where superscripts fc and fc — 1 refer to the different 
sides of the interface. The minus sign before the sec- 
ond term in the action reflects the opposite propagation 
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directions for the edge states of the two adjacent hq- 
uids. The operators exp(i<I>^) create equal charges for 
both values of p = /c, fc — 1. The prefactors in front of 
dt^^dx'^n ^'^^ opposite for different values of p. Hence, 
according to the criterion of stability of edge states^i^, 
the modes $^ and gap each other out for each 

n < k. Thus, the low-energy interface action reduces 
to L = l/(47r) / dtdx[K'^dt'Pldx^l-V{dx'^>l)^] and con- 
tains a single mode 



Let us now consider a system whose filling factor 
changes in a step- wise manner: as one moves across 
the edge of the quantum Hall bar, the filling factor 
first changes from to 1/3, then to 2/(2 -I- 2), then to 
3/(3 -|- 2) etc. The filling factor in the innermost part 
of the sample is k/{k + 2). The edge action can still 
be written in the form (|A3p . Different modes corre- 
spond to spatially separated interfaces between consec- 
utive regions with different filling factors. A quantum 
point contact in such systems brings close to each other 
two innermost edge channels, corresponding to the mode 

■ Clearly, any tunneling process can involve only those 
two interfaces. Hence, only one type of the quasiparti- 
cles with the Z-vector (0, ... ,0, 1) is allowed to tunnel. 
By shrinking the regions with intermediate filling factors 
< I' < k/{k + 2) one can deform the above system into 
a quantum Hall bar with a single interface between the 
1/ — and ly = k/{k + 2) regions. Still, it is clear that 
the innermost edge channel corresponds to the field 
The tunneling probability rapidly decreases as the path, 
traveled by the tunneling particle, grows. Hence, the 
maximal tunneling probability corresponds to quasipar- 
ticles of type (0, . . . , 0, 1) and it is reasonable to neglect 
all other tunneling processes. 



Certainly, this is an approximation. In general, the 
tunneling probabilities (|A2p depend on the parameters 
Arm and 5m- Note that one phase Sm can be ex- 
cluded by absorbing it into the Aharonov-Bohm phase 
due to the magnetic flux $. Two more parameters can 
be excluded by tuning gate voltages at the point contacts. 
Still 3{k — 1) fitting parameters remain. An expression 
with a large number of fitting parameters is of limited 
use. Fortunately, the model with only one type of tun- 
neling quasiparticles becomes exact, if the edges of the 
Mach-Zehnder interferometer correspond to the interface 
between v = k/{k + 2) and v = {k — l)/{k + 1) liquids. 
This can be accomplished if in addition to the gray (green 
online) region with the filling factor 2 + fc/(fc-|-2) in Fig. 1 
and the white region with the filling factor one creates 
a strip with the filling factor 2 + (fc — l)/(fc + l) along 
the edges. As shown in the next section, it is possible 
to derive an expression for the Fano factor without any 
fitting parameters in that case. This situation will be our 
main focus below. 



2. The case of only one flavor allowed to tunnel 



In this section we will omit the index m in Eq. (jA2[) 
since it can assume only one value. We will also set 
5m = since it can be absorbed in the Aharonov-Bohm 
phase. We will focus on the limit of low temperatures 
so that quasiparticles only tunnel from the edge with the 
higher potential to the edge with the lower potential. The 
tunneling probability depends only on the number n of 
the trapped quasiparticles and reads 



Pn — A[l + ucos( 



27r$ 



(fc + 2)$o 



■Kn{k + 1) 
fc + 2 



)], (A5) 



The number of different probabilities p„ depends on the 
parity of k. For odd fc, Pn = Pn+k+2 and hence we can de- 
fine n modulo k+2. For even fc, p„ = Pn+2(k+2) Pn+k+2 
and n should be defined modulo 2(fc-|-2). Hence, the num- 
ber of different values of p„ equals ri,„ax = (fc -|- 2)[3 -I- 
(— 1)'^+^]/2. In what follows we will see a considerable 
difference between even and odd k. 

It will be convenient to use the following parametriza- 
tion of Eq. 



Pn — cxp(— i(/)o) + cxp(m0)] [v exp(i0o) -I- exp(— infl)], 

(A6) 

where A[l- \/l - u^]/2, v = [l + Vl - u^]/u, 0o = 
27r$/[(fc + 2)<I>o] and = 7r(fc + l)/(fc + 2). AU exponents 
exp(i7i0) are rimaxth roots of 1. 



a. Current 

We now use Eq. (|A6[) for the calculation of the current. 
The current / = limt_).oo Nq/t, where g = e/(fc -I- 2) is a 
quasiparticle charge and N is the large number of tun- 
neling events over a long time period t. It is easy to find 
the average time i„ between the tunneling events which 
change the number of the trapped quasiparticles from 
n — 1 to n and from n to n + 1: tn ~ 1/pn- Taking into 
account that there are only n„iax different values of p„, 
one finds the average time of rimax consecutive tunneling 
events: 



i,„ax lim t/N 

t-i-oo 



E 

n=l 



1/P« 



Hence, 



(A7) 



(A8) 



We next evaluate the sum i = E ^/Pn- Similar sums 
were evaluated numerically in the main part of the article. 
For an arbitrary k we need an analytic expression. Two 
different expressions hold for even and odd k. We will 
only discuss the details for even k as the case of odd k 
can be considered in a similar way. 
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Using the paramctrization (|A6p one gets 
,-_ E„ \Qn{vexp{-i(f>o))\^ 



B|P(uexp(-i0o))|' 



(A9) 



where P{z) = Hiiz + exp(iZ6l)) and Q„(z) = P{z)/{z + 

P{z) can be found from the basic theorem of algebra. 
Indeed, the roots of that polynomial are known and equal 
— exp{il6). There is only one such polynomial of power 
"■max ~ 2{k + 2) with the coefficient 1 before z""-". 
This polynomial is P{z) = z"""" — 1. From this one 
gets Q„(z) = Er=r~' z\~iy+^ cxpi-ien[l + 1]). The 
sum IQnP can now be calculated by first perform- 
ing the summation over n and reduces to ^„ |QnP = 
?^max(^'^"°""' — — !)■ Finally, for even k 



1) 



B(w2 _ 1)[„2 



cos nmax</ 



(AlO) 



quadratic and arithmetic means also implies that e* > q. 
The first inequality sets different upper limits on the ef- 
fective charge e* for even and odd k. For odd fc, e* < e. 
For even k, e* < 2e. This difference agrees with the dif- 
ference of the magnetic field dependences of the current, 
discussed above. Both upper limits can be reached as we 
will see below. 

As in the calculation of the current, we focus on the 
case of even k and only give the final answer for odd 
k. We need to find St^. One can easily see that in the 
notation of the previous subsection 



6t 



2 _ Y.n=C^ ^ |g„(llCXp(-Z0o))|4 

B^\P{vexp{-i(j>oW 

' ' l,m,p,q=0 

(A14) 



^S{p + q - I - m ~ n„iaxs), 



eB(w2 - 1) v^''+^ - 2u2fc+4 cos ^ 1 



V 



4fc+8 



1 



For odd fc, a similar calculation yields 



eBiv^ - 1) + 2v 



„k+2 



COS 



$0 



V 



2k+4 



(All) 



(A12) 



The expressions look similar for even and odd k but 
there is a significant difference between them. Indeed, for 
odd fc, the current is a periodic function of the magnetic 
flux with the period $o- For even k the period is two 
times shorter. Such superconducting periodicity reflects 
Cooper pairing of composite fermions. 



b. Noise 

The noise is defined as = 21imt^oo(Q^ — Q'^)/t, 
where Q is the total charge transmitted during the time 
period t. Let us set t = mi, where m is a large integer. 
Such choice of t corresponds on average to N = mrimax 
tunneling events. The total time required for N tunnel- 
ing events can be expressed as r = mt + St, where St is 
a fluctuation. The total charge transferred through the 
interferometer after N events is exactly Nq. In a good ap- 
proximation, the charge transferred over the time t = mt 
is then Q = Nq — ISt, where / is the average current 
(|AlllAl2p . Substituting this expression in the deflnition 
of the noise one gets S — 2I^ST^/t. The calculation of 
St"^ is straightforward. One flnds: St^ = mSt^ , where 
St^ = Yl=7 ^Ivi- Finally, the Fano factor 



S_ 
21 



(A13) 



Some restrictions on e* are evident from elementary 
inequalities. Obviously, e* < rimax'?- The inequality of 



where s is an integer and the discrete delta function 
S(0) — 1, S{a 7^ 0) = 0. The remaining summation is 
tedious but straightforward. For even k we obtain 



" k + 2^ yo-2k+i -i' ^dw dw w -I 

(4fc + 7 - w—f -h 

aw w — 1 

2w^+^ cos(-— )(2A; + 2,- w — )- ],(A15) 

$0 dw dw w — 1 



where w ~ [1 + Vl — v?]'^ /v?. 

A similar calculation for odd k yields 



e U) - 1 2r ^ ^ w^^^ - 1 
k + 2^w^+'^ - V ^dw^dw w-1 



{2k + 3-w-^f 
dw 



w — 1 



2w^^^^>'^ cos{——){k + 1 - w—)- ](A16) 

<f>o dw dw w ^ \ 

The expressions are rather similar for even and odd k 
but their periodicity as a function of the magnetic flux is 
different just like for the current. 

For a general u, the above expressions are complicated. 
They simplify in an important limit case considered be- 
low. By tuning the gate voltages at the tunneling con- 
tacts it is always possible to make equal the tunneling 
amplitudes at the two contacts. The desired situation 
can be achieved by selecting such gate voltages that the 
total currents would be the same when only one contact 
is open no matter which one it is. As discussed in section 
V, this corresponds to u = 1 at a sufficiently low tem- 
perature and voltage bias. The expressions for the Fano 
factor greatly simplify in that limit. For even k, 



8(fc + 2)2 



1 + [4(fc + 2)2-1] cos ^ 



6(A: + 2) 



(A17) 
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For odd k, 



2{k + 2f + 1 - [(fc + 2)2 - 1] cos 



27r'I> 



3(fc + 2)2 



(A18) 



3 ^ 

and 



As a function of the flux, e* oscillates between e[ 

3(fc+2)2 ] ^'^'^ oscillates between + 3(^.^2)2 ] 

e. The maximal values of the Fano factor are thus 2 x e 
at odd k and 2 x 2e at even A:. The minimal values of the 
Fano factor uniquely identify states with different k. 



3. General case 

Here we address the situation when all quasiparticle 
flavors are allowed to tunnel. As discussed above, this 
situation is less interesting than the case of just one flavor 
allowed to tunnel. Indeed, the expressions for the current 
and noise depend on numerous fitting parameters and 
thus are less useful than the results (|A17IA18p . Besides, 
in a general case it would be difficult to tune the system 
to reach its maximal possible Fano factor. Instead, there 
is an interval of relevant Fano factors and this makes 
the identification of a state more difficult. Nevertheless, 
we discuss below how one can calculate the transport 
properties for a general k when all flavors are allowed to 
tunnel. 

The current can be found from the system of equa- 
tions of the form Eq. (|7I8|) . /; in these equations are the 
probabilities to find the trapped topological charge I in 
the interferometer. Thus, it is important to understand 
how many different values of the topological charge arc 
possible. 

Different topological charges correspond to different rip 
in tunneling probabilities (jA2[) . Different rip may how- 
ever describe the same topological state. This happens, if 
all tunneling probabilities (|A2[) are the same for a certain 
set of Up and another set of Up -\- Ap. The probabilities 
are equal, if for each m 



7tA,„ , = 27rr„ 



(A19) 



where r„j is an integer. One finds from (|A19|1 that A„i = 
2?^™, + J2p ^pl (fc + 2). Adding up such equations for aU 



TO, one finds that Ap = (fc -|- 2) rp and hence 



(Ai,...,Afc)= ^ 



(A20) 



where d„j = (1, . . . , 1, 3, 1, . . . , 1) and the number 3 stays 
in position to. Eq. (jA20|l means that adding or subtract- 
ing a vector from the ^-vector n = (ni, . . . ,nk) does 
not change the topological charge. d„i can be understood 
as Z- vectors of electrons. 

If wc choose ri^ — 1 and rp = — 1 for one p < k, then 
adding the ^-vector (jA20p to n increases rij, + 2 

and decreases rip Up ~ 2. Setting = — 1 and rp = 1 



corresponds to the operation ^ — 2, rip rip 2. 
This operation involves adding to n one vector d^ and 
subtracting another vector d„. Performing several such 
operations, one can always reduce all rip with p < fc to 
zeroes and ones. Choosing = fc -I- 1 and rp = — 1 for all 
p < k corresponds to the operation — > rik + 2(fc + 2), 
rip — Hp. The operation involves 2(fc + 1) additions and 
subtractions of d^'s. Setting r^ = — (fc + 1) and Vp = 1 
for all p < fc corresponds to —> — 2 (fc + 2), rip —> 
rip. Hence, by adding and subtracting vectors d^ one 
can always reduce Z-vectors to topologically equivalent 
vectors of the form (si, . . . , Sk), where Sp = 0, 1 at p < fc 
andO < Sk < 2(k+2). We will call such vectors allowable. 
The total number of the allowable vectors is 2'^(fc + 2). 

We have established that any /-vector can be reduced 
to the allowable form by an even number of additions and 
subtractions of vectors d^. Moreover, one can get ex- 
actly one allowable vector Ne(n) from each n by an even 
number of additions and subtractions. Indeed, an even 
number of additions and subtractions does not change 
the parity of each component of the /-vector. This fixes 
Sp, p < k. The residue rim) mod 2(fc-|-2) is also invari- 
ant with respect to an even number of subtractions and 
additions. This fixes Sk- As a consequence, no other al- 
lowable vectors can be obtained from an allowable vector 
by an even number of additions and subtractions. 

Next, note that No(n) = Ne(n -I- di) is the only al- 
lowable vector that can be obtained from n by an odd 
number of additions and subtractions. Indeed, if another 
allowable vector Mo(n) could be obtained from n by an 
odd number of additions and subtractions then No(n) 
and Mo(n) could be obtained from each other by an even 
number of additions and subtractions: one first gets n 
from Mo(n) and then No(n) from n. We have already 
proved that this is impossible. 

Thus, exactly one allowable vector Ne(n) can be ob- 
tained from n by an even number of additions and sub- 
tractions and exactly one allowable vector No(n) can be 
obtained by an odd number of additions and subtrac- 
tions. Note that Ne(n) 7^ No(n) since the parity of each 
component Sp changes after each addition or subtraction 
and hence the parities of the components of Ne(n) and 
No(n) arc opposite. Note also that Ne(n) = No(No(n)) 
and No(n) = No(Ne(n)). Thus, the set of allowable vec- 
tors consists of 2*^"^ (fc -1-2) pairs of the form (V,No(V)). 
Different pairs have no common elements, /-vectors from 
different pairs cannot be obtained from each other by any 
number of additions and subtractions. If an /-vector can 
be reduced to the allowable vector V by an even num- 
ber of additions and subtractions it can also be reduced 
to No(V) by an odd number of additions and subtrac- 
tions. Let us now select one arbitrary vector from each 
pair of allowable vectors (V,No(V)). We get a set H 
of 2''"~-'^(fc + 2) allowable vectors. An arbitrary /-vector 
(rii, . . . , rik) is topologically equivalent to one of the vec- 
tors in J. Hence, there are 2'^"^(fc + 2) different topolog- 
ical charges. 

Introducing a probability distribution /; for different 
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trapped topological charges /, one can compute the cur- 
rent from the system of equations (|7|8p . Even for fc = 2 
the result is complicated and in the main part of the ar- 
ticle we showed it only in the graphical form (Figs. 3, 
4) except for the flavor-symmetric case when it simplifies 
considerably. 

A simplification is also possible in the flavor-symmetric 
case for a general k. In that case one can remove 
the index m from the constants Am, Um and 6m since 
they do not depend on the flavor. Let us set n ~ 
{J2 "p) 2(A:-|-2) and S = X]("-p ^^od 2), i.e., let n de- 
note the total number of trapped quasiparticles modulo 
2{k + 2) and S show how many rip's are odd. Note that 
n and S have the same parity. The number of possible 
pairs (?!, S) equals C = (fc-|-l)(fc-|-2). Let us introduce a 
distribution function fn,s- A simplification of the steady 
state equations results from the fact that one can write 
a closed set of equations for fn,s- 



n,S 



dt 

-A[l + u cos( 



Ml- 



u cos( 



27r$ 



where the notation implies that fn,s = for S* < and 
S > k. The current can then be expressed as 



fn,s{A[l - ucos( 



27r<I> 



(fc 



+A[l + ucos( 



27r$ 



[k + 2)$o 



+ 5- 



5- 



fc + 2 



)](fc-5m22) 



System (jA2ip contains C = (fc + l)(fc -I- 2) equations, 
much fewer than (fc + 2)2*^^^ in a general case. More- 
over, one can reduce the total number of equations to 
no more than k/2+ 1. Indeed, Eq. (jA21[) allows one to 
express fn,s via /„-i.s±i- /,i-i,s±i can be expressed via 
/„_2,s', where S" = - 2, S*, + 2, etc. After 2(fc + 2) 
steps, one expresses fn.s via the values of the distribu- 



(fc 



-2)$o 
7rn 



fc + 2 



)]'$ion function /„_ 



where 5" = 0, . . . , fc. Since n is 



-l-/„_i,s-i(fc - S - 



(fc + 



2)$o 
I- u cos( 



fc4 
27r$ 



l)A[l - MC0S( 



(fc + 2)$o 
(fc + 2)$c 



+ 5- 



+ 5- 



-2(k+2),Si 

defined modulo 2(fc -I- 2), this means that a closed system 
S)} X /„ ^an be obtained for fc + 1 variables fn,s for any fixed n. 

' At least fc/2 of those variables are zeroes since n and S 
7r(n — 1) have the same parity for any nonzero fn,s- This leaves no 
'^niore than fc/2 -I- 1 equations. After the system is solved, 
■nin — 1] ^^"^ '^^'^ immediately compute fn+i,s from /„ 5 with Eq. 
f )| <pi2ip . then express /„+2,s via /„+i,s and so on. 



^ The Quantum Hall Effect, edited by R. E. Prange and S. M. 

Girvin (Springer- Verlag, New York, 1987). 
^ C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. 

Das Sarma, Rev. Mod. Phys. 80, 1803 (2008). 
^ A.Yu.Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 
* For a review of proposed states see B. J. Overbosch and 

X.-G. Wen, arXiv:0804.2087 (unpublished). 
^ G. Moore and N. Reed, Nucl. Phys. B 360, 362 (1991). 
® S.-S. Lee, S. Ryu, C. Nayak, and M. P. A. Fisher, Phys. 

Rev. Lett. 99 236807 (2007). 

M. Levin, B. I. Halperin, and B. Rosenow, Pliys. Rev. Lett. 

99, 236806 (2007). 
® B. I. Halperin, Helv. Phys. Acta 56, 75 (1983). 
^ X.-G. Wen and A. Zee, Phys. Rev. B 46, 2290 (1992). 
1° J. Yang, Z. Su, and W. Su, Mod. Phys. Lett. B 6, 119 

(1992). 

" B. Blok and X.-G. Wen, Phys. Rev. B 42, 8145 (1990). 
^2 R. H. Morf, Phys. Rev. Lett. 80, 1505 (1998). 

E. H. Rezayi and F. D. M. Haldane, Phys. Rev. Lett. 84, 

4685 (2000). 

" G. Moller, S. H. Simon, Phys. Rev. B 77, 075319 (2008). 
A. E. Feiguin, E. Rezayi, K. Yang, C. Nayak, and S. Das 
Sarma, Phys. Rev. B 79, 115322 (2009). 
C. Toke, N. Regnault, and J. K. Jain, Phys. Rev. Lett. 98, 
036806 (2007). 

" G. R. Dean, B. A. Plot, P. Hayden, S. Das Sarma, G. 

Gervais, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 

101, 186806 (2008). 
1* G. Zhang, T. Knuuttila, Y. Dai, R. R. Du, L. N. Pfeiffer, 



and K. W. West, Phys. Rev. Lett. 104, 166801 (2010). 
M. Stern, P. Plochocka, V. Umansky, D. K. Maude, M. 
Potemski, and L Bar- Joseph, Phys. Rev. Lett. 105, 096801 
(2010). 

^° A. Wojs, G. Moller, S. H. Simon, and N. R. Cooper, Phys. 

Rev. Lett. 104, 086801 (2010). 

S. Das Sarma, G. Gervais, and X. Zhou. ' arXi^:1007.1688l 

(unpublished) . 
22 J. K. Jain, Physics 3, 71 (2010). 

M. R. Peterson and S. Das Sarma, Phys. Rev. B. 81, 

165304 (2010). 

I. P. Radu, J. B. Miller, C. M. Marcus, M. A. Kastner, 
L. N. Pfeiffer, and K. W. West, Science 320, 899 (2008). 
A. Stern, Nature 464, 187 (2010). 

E. Fradkin, C. Nayak, A. Tsvelik, and F. Wilczek, Nucl. 
Phys. B 516, 704 (1998). 

S. Das Sarma, M. Freedman, and G. Nayak, Phys. Rev. 
Lett. 94, 166802 (2005). 

A. Stern and B. L Halperin, Phys. Rev. Lett. 96 016802 
(2006). 

P. Bonderson, A. Kitaev, and K. Shtengel, Phys. Rev. Lett. 
96, 016803 (2006). 
^° R. L. Willett, L. N. Pfeiffer, and K. W. West, Proc. Natl. 
Acad. Sci. U.S.A. 106, 8853 (2009). 

P. Bonderson, C. Nayak, and K. Shtengel, Phys. Rev. B 
81, 165308 (2010). 

A. Stern, B. Rosenow, R. Ilan, and B. Halperin, 

,iarXnFy09 12.4394 (unpublished). 

R. Ilan, B. Rosenow, and A. Stern, larXiv: 1005.47721 (un- 



13 



published). 

Y. Ji, Y. C. Chung, D. Sprinzak, M. Heiblum, D. Mahalu, 
and H. Shtrikman, Nature 422, 415 (2003). 
K. T. Law, D. E. Feldman, and Y. Gefen, Phys. Rev. B. 
74, 045319 (2006). 

D. E. Feldman and A. Kitaev, Phys. Rev. Lett. 97, 186803 
(2006). 

D. E. Feldman, Y. Gefen, A. Kitaev, K. T. Law, and A. 

Stern, Phys. Rev. B 76, 085333 (2007). 

K. T. Law, Phys. Rev. B 77, 205310 (2008). 

V. V. Ponomarenko and D. V. Averin, Phys. Rev. B 80, 

201313 (2009). 

X. -G. Wen, Quantum Field Theory of Many-Body Systems 
(Oxford University Press, Oxford, 2004). 



M. Dolev, M. Heiblum, V. Umansky, A. Stern, and D. 
Mahalu, Nature (London) 452, 829 (2008). 
A. Kitaev, Ann. Phys. (N.Y.) 321, 2 (2006). 
N. Byers and C. N. Yang, Phys. Rev. Lett. 7, 46 (1961). 
R. de Picciotto, M. Reznikov, M. Heiblum, V. Umansky, G. 
Bunin, and D. Mahalu, Nature (London) 389, 162 (1997). 
L. Saminadayar, D. G. Glattli, Y. Jin, and B. Etienne, 
Phys. Rev. Lett. 79, 2526 (1997). 

E. B. Fel'dm an, Theor. Exp. Chem. 10, 645 (1976). 
M. Heiblum. la"rXiv:091 2.4868 (unpublished). 

F. D. M. Haldane, Phys. Rev. Lett. 74, 2090 (1995). 
H.-c. Kao, G.-H. Ghang, and X.-G. Wen, Phys. Rev. Lett. 
83, 5563 (1999). 



